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Narrowband Angular Reflectance Properties of 
the Alkali Flats at White Sands, New Mexico 


Charles H. Whitlock * Staart R. LeCroy , f and Robert J. Wheeler f 


estilts from helicopter measurements of the angular 
properties of surface reflectance for the alkali flats regions 
of the White Sands Missile Range are presented for the 
wavelength interval of 0.4 to 0.85 pm. This work was 
performed to allow accurate radiative transfer calcula- 
tions over the region. Detailed tables and interpolation 
equations are given that permit other investigators to 
perform satellite calibrations over the alkali flats site. The 
effects of wavelength and soil moisture on narrowband 
angular reflectance are also investigated. Although there 
is a spectral variation in surface albedo, there is little 
spectral effect in Anisotropic Factor except in the forward- 
scattering peak at solar zenith angles greater than 60°. 
The magnitude of the forward -scattering peak is also 
sensitive to soil moisture, with wet conditions causing a 
larger peak. The significance of this result is that angular 
reflectance properties at the center of the alkali flats 
usually will be different than those at the flats edge 
because moisture differences typically exist. 


INTRODUCTION 

Although satellite instruments are usually calibrated 
before flight, there is generally a need for recalibration 
once flight conditions have been achieved. These calibra- 
tions are often done using homogeneous, bright ground 
targets such as deserts. The ground-target approach 
(indirect method) employs a radiative transfer model 
to determine the upwelling directional radiance that 
reaches the satellite instrument’s detector (Slater et 
ah, 1987; Frouin and Gautier, 1987). This calculated 
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radiance is then compared to the digital instrument 
counts to determine the calibration coefficients. The 
success of this technique depends on how accurately 
both surface reflectance properties and atmospheric 
characteristics can he characterized. For bright targets, 
one key parameter is the angular distribution of surface 
reflectance (Jackson ct ah, 1990). This paper will only 
consider the surface reflectance as a function of solar- 
zenith, view-zenith, and view-azimuth angles for the 
wavelength interval 0.4 to 0.85 /mi. 

Early satellite calibration studies (Frouin and Gau- 
tier, 1987) used laboratory-derived angular reflection 
properties, based on samples of White Sands Missile 
Range (WSMR) gypsum sand (Walraven and Coulson, 
1972). Recent studies have shown that the angular re- 
flectance properties of the WSMR dunes region are not 
the same as those of the WSMR alkali flats (Markham 
et ah, 1990; Deering et ah, 1990). Broadband (0. 2-4.0 
gm) angular reflectance measurements of the alkali flats 
region have been made (Eaton and Dirmhirm, 1979), 
but narrowband measurements at the edge of the flats 
indicated the possibility of both spectral and surface 
texture effects (Deering et ah, 1990). 

From this history, it was deemed necessary to per- 
form new narrowband angular reflectance measurements 
over the central part of the WSMR alkali flats region. 
This region was selected because it was a large, low- 
vegetated, semi-homogenous area suitable for use as a 
calibration target for operational satellites with pixel 
sizes up to 3 km. This paper presents dedicated mea- 
surements that will permit other investigators to per- 
form satellite calibrations over the alkali flats site. 

Measurement Technique 

Measurements of angular reflectance variability have 
been accomplished from a number of platforms, such 
as small towers or booms (Eaton and Dirmhirm, 1979; 
Deering and Leon, 1986; Jackson et ah, 1990), fixed- 
wing aircraft (Davis and Cox, 1982; Markham et ah. 
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Figure 1. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.30 for a solar- 
zenith angle of 15°. 


1990; Irons et ah, 1991), and helicopters (Whitlock et 
ah, 1987; Williams et ah, 1991; Purgold et ah, 1994). 
The particular technique used depends on the target 
size, texture, and homogeneity combination; wavelength 
requirements; angular requirements; and level of sophis- 
tication of the radiative transfer model being used. 

In our case, the initial requirement was for angular 
properties of surface reflectance for a 3 x 3 km target 
near the center of the alkali flats at WSMR. The ultimate 
use of the data was in the Finite Difference Method 
radiative transfer model (Barkstrom, 1976; Suttles, 1981) 
to validate Advanced Very High Resolution Radiometer 
(AVHRR) and Visible Spin Scan Radiometer (VISSR) cali- 
bration coefficients. This model was selected for use 
because it considered all orders of scattering, has been 
intercompared with other radiative transfer methods (Sut- 
tles, 1985; Ellingson and Fouquart, 1990), and was verified 
with experimental data (Whitlock et al., 1985; 1990). The 
following angular reflectance cases were desired: 

1. Solar zenith angles from 15° to 80°. 

2. View-zenith angles from 0° to 75°. 


ZENITH 



Figure 2. Bidirectional characteristics of WSMK at a wave- 
length of 0.65 micron and an albedo oi 0.30 for a solar- 
zenith angle of 40°. 

3. View-azimuth angles from 0° to 180° relative to 
the plane of the sun. Symmetry is assumed about 
the plane of the sun. 

4. Selected spectral bands between 0.40 and 0.85 
jum with bandwidths less than 0.07 fim. 

5. A range of soil moisture conditions. 

The scanning-radiometer helicopter system described 
in Purgold et al. (1994) was used to perform the basic 
radiance measurements during May 1991. Simultaneous 
reflectance measurements were made at 0.4, 0.55, 0.65, 
and 0.85 gm. Bandwidths were approximately 0.07 gm. 
The partial-clovorleaf flight pattern required 8 minutes 
to measure surface radiance from 0° to 75° in view- 
zenith angle at eight different view-azimuth angles rela- 
tive to the solar plane. Patterns were flown for 33 
different solar-zenith angles and two different soil mois- 
tures. The helicopter radiance data were integrated over 
angle to determine surface albedo, validated against 
surface albedo measurements (Wheeler et al., 1994), 
and converted to Anisotropic Factors (AF) following 
Suttles et al. (1988). AF is defined as the ratio of fl 
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Figure 3. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.30 for a solar- 
zenith angle of 50°. 

times the directional radiance at a given view angle 
divided by the upward hemispherical irradiance from 
the surface. AF equals 1.0 at all angles for a Lambertian 
surface. Non-Lambertian AF values may be converted 
to relative bidirectional reflectance factor (RBRF) values 
(Jackson et ah, 1990) by dividing the AF value at an 
arbitrary view-zenith and view-azimuth angle by the AF 
value at nadir (view-zenith = 0). AF values are validated 
following a concept analogous to that of Moran et al. 
(1990). 

RESULTS 

The relation between soil moisture and surface albedo 
is unknown in this region because standard moisture 
test procedures do not provide consistent results with 
the alkali flats material. Standard bake-out temperatures 
do not remove all water in the gypsum / salt / algae / 
lichen soil mixtures. For this reason, we use surface 
albedo as an indicator of soil moisture variability in this 
report. Wheeler et ah (1994) indicate that surface albedo 
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Figure 4. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.30 for a solar- 
zenith angle of 60°. 


values range from 0.22 (saturated with water) to 0.52 
(unusually dry conditions) at a wavelength of 0.55 gm. 
Typical albedo values during the months of January, 
April, May, November, and December (when minimum 
cloud cover is observed) range from 0.30 to 0.45 with 
the variation on a daily basis depending on local rainfall. 
WSMR dune's albedos are higher but with a similar 
range of variability. 

All surface albedo values are given for a wavelength 
of 0.55 pm for the remainder of this report. This wave- 
length is used because inexpensive instrumentation is 
commercially available (LI-COR Model 210SB) for mon- 
itoring albedo from 0.5 to 0.6 p m. Values at 0.40, 0.65, 
and 0.85 pm wavelengths may be obtained by multi- 
plying the 0.55 pm value by the following correction 
factors: 

Wavelength, A 0.40 pm 0.55 pm 0.65 pm 0.85 pm 

Correction Factor 0.48 1.00 1.22 1.26 

The above correction factors are the average of a num- 
ber of albedo measurements at each of the four wave- 
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Figure 5. Bidirectional characteristics of VVSMH at a wave- 
length of 0.65 micron and an albedo of 0.30 for a solar- 
zenith angle of 70°. 


lengths at two albedo values from the helicopter. Mea- 
surements at albedos of 0.30 and 0.42 indicate that 
moderate variations in soil moisture do not have a 
significant effect on spectral characteristics of surface 
albedo for the alkali flats. 

For reasons of brevity, results have been synthe- 
sized to a common set of solar-zenith, view-zenith, and 
view-azimuth angles where straight-line interpolations 
may be used to represent the data. Figures 1 to 6 show 
AF results at 0.65 fim for a surface albedo of 0.30 (at 
0.55 fim) for six solar zenith angles. These data represent 
moderately damp soil-moisture conditions. Figures 7 to 
12 show the same results for a surface albedo of 0.42 
representing moderately dry conditions. Data for 0.65 
fim were selected for display because that wavelength 
is most common to many different satellite instruments. 

AF values for other wavelengths within the range 
of 0.40 to 0.85 fim may be computed from: 

AF; = AFa= ( i.65*[1 + (A - 0.65)*TAN p\ (1) 
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Figure 6. Bidirectional characteristics of WSMK at a wave- 
length of 0.65 micron and an albedo of 0.30 for a solar- 
zenith angle of 80°. 

w'here: 

A = wavelength of interest in fim. 

AF ; =<>.65 = value from Figures 1 to 12 for the desired 
solar-zenith, view-zenith, view-azimuth, 
and surface albedo combination. 

P = value from Table 1 for the desired 

solar-zenith, view-zenith, view-azimuth, 
and surface albedo combination, p is the 
slope of a least-squares-fit line on a plot 
of the AF; to AFA = 0.65 pratio as a 
function of A. 

In general, the relative accuracy of the calculated AF; 
is within 6% of the measured AF; for all solar-zenith, 
view -zenith, view'-azimuth, and surface albedo combina- 
tions except portions of the forward-scattered peaks 
at the 70° and 80° solar-zenith angles. The following 
correction to equation (1) is required for those view- 
zenith and view-azimuth angles within the forward- 
scattered peak: 
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Figure 7. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.42 for a solar- 
zenith angle of 15°. 


AFa = £*[AFa,eqn(1)] (2) 

where £ is a correction factor given in Tables 2 and 3 
for the angles within the forward-scattering peak. 

It is necessary to adjust AF values to account for 
the variation in soil moisture. The authors use a linear 
interpolation / extrapolation process based on surface 
albedo values at 0.55 fim. The linear process is necessary 
because angular measurements were made at only two 
values of surface albedo. The adjustment to obtain AF 
values at surface albedo, A, is: 

AF a = AFa-0.30 + [(A — 0.30)*(AFo.42 — AFo 30 ) / 0.12] (3) 

As the soil becomes more saturated with water, specular 
reflections effects may alter AF values in a nonlinear 
manner. Therefore, the authors do not recommend us- 
ing the AF data in this report for surface albedos below 
0.27. 
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Figure 8. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.42 for a solar- 
zenith angle of 40°. 


Validation 

The above helicopter-derived results may be verified in 
a limited manner by using satellite data from different 
view angles similar to Moran et al. (1990). Satellite 
digital counts should be approximately proportional to 
the radiance from the surface in the direction of the 
satellite if the atmospheric variation is not large and 
satellite calibration does not vary by a large amount. 
Directional radiance at the surface is proportional to 
scaled directional reflectance, which is the surface al- 
bedo times AF scaled by the cosine of the solar-zenith 
angle (SZA). Figure 13 shows NOAA-11 band-1 counts 
as a function of surface scaled directional reflectance 
calculated from surface-site albedo data (Wheeler et al., 
1994) and helicopter AF values for a 31-month period 
beginning November 1988. 

Satellite counts varied by a factor of two, depending 
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Figure 9. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.42 for a solar- 
zenith angle of 50°. 


Figure 10. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.42 for a solar- 
zenith angle of 60°. 


mostly on a wide range of combinations of surface 
albedo, solar zenith angle, and AF, which existed during 
the period. Surface albedo ranged from 0.27 to 0.51, 
and the solar-zenith angle varied from 18° to 61°. 
View-zenith and view-azimuth angles ranged from 1° 
to 41° and 2° to 171°, respectively. The solar-zenith, 
view-zenith, and view-azimuth combinations caused the 
AF to range from 0.89 to 1.07, with a 20% variation. 

Aerosol optical depth varied between 0.012 and 
0.132 (LeCroy et ah, 1994), and water vapor ranged 
from 0.1 to 0.8 cm from surface-site data. Radiative 
transfer calculations with these inputs (assuming Conti- 
nental aerosols) indicate only 5% change in top-of- 
atmosphere radiance (or satellite counts). This uncer- 
tainty would show up as random scatter about the 
regression line in Figure 13. 

Che and Price (1992) suggest that AVHRR band 1 
may have slowly degraded as much as 3.5% per year 
(with absolute uncertainty equal to 7%). More recent 
results (Rao and Chen, 1993) indicate a 1.2% per year 


degradation rate. Analysis of the 32 cases plotted in 
Figure 13 does not indicate a large degradation trend. 
Points separated by a 30-month time difference fall 
within 3% of each other for similar values of scaled 
directional reflectance. Random scatter caused by aero- 
sol and water vapor variation may be partially hiding 
calibration degradation, however. 

The scatter of the maximum counts about the re- 
gression line in Figure 13 is 6.1% (rms error equals 
3.3%). Because this is similar to the atmospheric and 
calibration uncertainty, measurement errors for both 
surface albedo and AF are probably less than 6% . There 
is no indication that either high or low values of AF 
cause counts to be either above or below the line for 
the 32 overflights in this experiment. Likewise, there is 
no evidence that the linear interpolation / extrapolation 
process for soil moisture effects [Equation (3)] is not 
adequate for the albedo range of these overflights. Pairs 
of points at both ends of the albedo range (0.27 to 0.51) 
fall very close to the line at similar values of scaled 
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Figure 11. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.42 for a solar- 
zenith angle of 70°. 


Figure 12. Bidirectional characteristics of WSMR at a wave- 
length of 0.65 micron and an albedo of 0.42 for a solar- 
zenith angle of 70°. 


directional reflectance. Results from this analysis give 
increased confidence in both the helicopter-measured 
AF values and the linear correction for soil moisture 
effects. 

It should be noted that a significant change in 
aerosol type or magnitude can alter the validity of the 
above test. The data shown in Figure 3 are for the usual 
aerosol conditions at WSMR prior to arrival of the Mt. 
Pinatubo volcanic plume. Satellite counts taken within 
6 months after the volcanic cloud reached WSMR were 
as much as 15% lower than would be expected, and 
the slope of the regression line was different than that 
in Figure 13. At some times during the post-Pinatubo 
period, volcanic-cloud optical depths in the stratosphere 
w ere five times that typical of the lower-altitude aerosol 
column. Knowledge of aerosol optical depth magnitudes 
and aerosol type is required to assess uncertainties 
caused by aerosols in the analysis. 


ANALYSIS 

It is of interest to assess both spectral and soil moisture 
effects on the angular reflectance properties of the alkali 
flats. For ease of comparison, we convert AF into RBRF 
in the plane of the Sun using the procedure discussed 
previously. Figures 14 and 15 show RBRF values for a 
SZA of 60° at the tw ? o albedo values for which helicopter 
measurements were made. In both cases, the major 
wavelength effect appears to be in the peak in the 
forward-scattering direction. The cause of this charac- 
teristic has not been established. Table 1 shows low 
values of fl at low values of SZA indicating that wave- 
length effects become smaller as the AF forward-scatter- 
ing peak is reduced in size. 

Soil moisture effects are examined at two locations 
on the alkali flats region. As noted earlier, the helicopter 
measurements were conducted over a 3x3 km central 
flats test area (CFTA) in May 1991. On November 16, 




NOAA-1 1 COUNTS (1 0-BIT) 


178 Whitlock et al. 


Table 1. Beta Angles for Anisotropic Factor Spectral Correction 


Solar 

zenith 

View 

zenith 

View azimuth (0-45) 

View 

azimuth 

(90) 

View 

azimuth 

(135-180) 

Albedo -0.42 

Albedo = 0.30 

Both Albedos 

15 

0 

1 

2 

2 

2 


15 

1 

1 

1 

2 


30 

1 

0 

1 

2 


45 

1 

- 1 

0 

I 


60 

0 

-3 

- 1 

0 


75 

- 1 

-6 

- 1 

- 1 

40 

0 

3 

6 

4 

5 


15 

2 

3 

4 

5 


30 

1 

0 

3 

4 


45 

- 1 

-3 

1 

2 


60 

-2 

-9 

- 1 

- 1 


75 

-4 

- 15 

-3 

-3 

50 

0 

3 

7 

5 

6 


15 

3 

4 

5 

6 


30 

2 

0 

4 

5 


45 

-2 

-3 

3 

3 


60 

-6 

-10 

1 

0 


75 

- 11 

- 18 

-5 

- 1 

60 

0 

5 

5 

5 

5 


15 

0 

5 

4 

6 


30 

-5 

2 

3 

7 


45 

-9 

-5 

1 

6 


60 

-12 

- 19 

0 

5 


75 

- 15 

-21 

- 1 

4 

70 

0 

4 

10 

6 

10 


15 

-4 

5 

7 

10 


30 

- 16 

- 1 

6 

11 


45 

- 17 

-17 

4 

12 


60 

-16 

-32 

2 

11 


75 

- 12 

-39 

- 1 

10 

80 

0 

6 

10 

9 

9 


15 

-6 

2 

8 

13 


30 

- 18 

-7 

7 

16 


45 

-31 

-31 

6 

17 


60 

-35 

-48 

3 

19 


75 

-30 

- 47 

2 

20 

All values 

» are in degi 

rees, and 0° azimuth 

is toward the Sun 

in the solar pi; 

ane. 


Figure 14. Relative Bidirectional Reflectance Function for 
albedo of 0.30 and solar-zenith angle of 60°. 


Figure 13. NOAA-1 1 band-1 counts versus scaled direc- 
tional reflectance of surface. 
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Table 2. Epsilon Correction Factor for Forward-Scattered Peak for Surface 
Albedo = 0.30 


Solar 

zenith 

View 

azimuth 

View 

zenith 

A = 0.4 ju?n 

X = 0. 55 nm 

A = 0.85 fwi 

70° 

0° 

30° 

0.98 

1.00 

0.99 



45° 

0.94 

0.99 

0.95 



60° 

0.90 

1.00 

0.91 



75° 

0.83 

0.98 

0.93 

70° 

25° 

30° 

0.99 

0.99 

0.99 



45° 

1.03 

0.99 

0.99 



60° 

0.98 

0.93 

0.92 



75° 

0.91 

0.95 

0.84 

70° 

45° 

30° 

1 .00 

1.00 

1.00 



45° 

1.02 

0.99 

1.01 



60° 

0.97 

0.93 

0.96 



75° 

0.92 

0.95 

0.89 

80° 

0° 

30° 

0.93 

0.96 

0.92 



45° 

0.88 

0.92 

0.87 



60° 

0.80 

0.89 

0.84 



75° 

0.75 

0.85 

0.81 

80° 

25° 

30° 

0.93 

0.95 

0.91 



45° 

0.93 

0.92 

0.90 



60° 

0.84 

0.84 

0.84 



75° 

0.78 

0.83 

0.74 

80° 

45° 

30° 

0.94 

0.96 

0.95 



45° 

0.92 

0.92 

0.94 



60° 

0.85 

0.85 

0.89 



75° 

0.80 

0.83 

0.78 


Table 3. Epsilon Correction Factor for Forward-Scattered Peak for Surface 
Albedo = 0.42 

Solar 

zenith 

View 

azimuth 

View 

zenith 

A - 0.4 fim 

A = 0.55 |/wi 

A = 0.85 fim 

70° 

0° 

30° 

1.03 

0.97 

1.01 



45° 

1.01 

0.93 

0.97 



60° 

1.00 

0.91 

0.95 



75° 

0.99 

0.89 

0.93 

70° 

25° 

30° 

1.03 

0.97 

1.01 



45° 

1.01 

0.93 

0.97 



60° 

1 .00 

0.91 

0.95 



75° 

0.99 

0.89 

0.93 

70° 

45° 

30° 

1.01 

0.98 

1.02 



45° 

1.01 

0.94 

0.99 



60° 

1.00 

0.92 

0.97 



75° 

1.00 

0.90 

0.94 

80° 

0° 

30° 

0.95 

0.93 

0.92 



45° 

0.92 

0.88 

0.87 



60° 

0.87 

0.83 

0.82 



75° 

0.83 

0.77 

0.76 

80° 

25° 

30° 

0.95 

0.93 

0.92 



45° 

0.92 

0.88 

0.87 



60° 

0.87 

0.82 

0.82 



75° 

0.83 

0.77 

0.76 

80° 

45° 

30° 

0.95 

0.94 

0.96 



45° 

0.92 

0.89 

0.91 



60° 

0.88 

0.84 

0.85 



75° 

0.84 

0.79 

0.79 


80' 
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ALBEDO=0.42 
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Figure 15. Relative Bidirectional Reflectance Function for 
albedo of 0.42 and solar-zenith angle of 60°. 
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Figure 17 . Relative Bidirectional Reflectance Function for 
two locations on the alkali flats. Albedo values are for wave- 
length of 0.55 ,um, 


1988, Deering et al. (1990) performed similar measure- 
ments in a smaller area on the eastern edge of the flats 
at a location known as Chuck Site. Figure 16 illustrates 
the location of both areas. The flats at Chuck Site are 
at a higher altitude (5 to 15 m) and do not retain 
moisture as long as the CFTA. Chuck Site soil does not 
appear to contain as much algae and lichen as the 
CFTA, but surface texture is similar for both areas. 
Figure 17 compares RBRF for the two areas for several 
different albedo values. For the CFTA, it appears that 
the magnitude of the forward-scattering peak is a strong 
function of soil moisture as indicated by albedo value. 
More moist conditions with lower albedos have larger 
forward-scattering peaks. By coincidence, the Chuck 
Site albedo taken at the time of the Deering measure- 
ments in 1988 is close to the dry CFTA value in 1991. 
The helicopter CFTA measurements are in close agree- 
ment w ith the Chuck Site data of Deering et al. (1990). 
Precise agreement is not to be expected because of 
different measurement techniques and possible small 


Figure 16. White Sands /Tularosa Basin (Not to Scale). 



San Andres 
Mountains 


differences in surface texture and the amount of thin 
unstabilized gypsum dust (Wheeler et ah, 1994). 

CONCLUDING REMARKS 

Soil moisture causes significant change in surface albedo 
on the alkali flats region of WSMR. Increased wetness 
causes reduced surface albedo values at all wavelengths. 
The spectral variation of albedo relative to that at a given 
wavelength is lightly influenced by moderate changes 
in soil moisture, however. This suggests that albedo 
measurements at one wavelength may be used to esti- 
mate albedo values at other wavelengths within the 
interval 0.40 to 0.85 jum for a moderate range of soil 
moisture variation. 

The anisotropic factor exhibits minimal spectral varia- 
tion except in the forward- scattering peak at significant 
solar-zenith angles. The magnitude of the forward-scat- 
tering peak is also sensitive to soil moisture, with wet 
conditions causing a larger peak. 

Comparison of helicopter-derived data with that of 
Deering et al. (1990) indicates that the two methods 
give nearly identical results for similar combinations of 
solar-zenith angle and surface albedo at two different 
locations on the alkali flats. 

Caution is suggested when comparing satellite counts 
of both the central flats and Chuck Site regions on the 
same satellite image for large solar-zenith and view- 
zenith angles. Typically, the central flats region has a 
low'er albedo than Chuck Site because of lower elevation 
and higher soil moisture. Under such conditions, mois- 
ture variation will cause the anisotropic factor to be 
significantly different for each of the two regions. Satel- 
lite radiance differences will no longer be directly pro- 
portional to surface albedo differences, as is sometimes 
assumed. This issue is minimized for low solar-zenith 
and view-zenith angle conditions. 
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